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ABSTRACT 

Sodium abundances have been determined in a large number of giants of open clusters 
but conflicting results, ranging from solar values to overabundances of up to five 
orders of magnitude, have been found. The reasons for this disagreement are not well- 
understood. As these Na overabundances can be the result of deep mixing, their proper 
understanding has consequences for models of stellar evolution. As discussed in the 
literature, part of this disagreement comes from the adoption of different corrections 
for non-LTE effects and from the use of different atomic data for the same set of lines. 
However, a clear picture of the Na behaviour in giants is still missing. To contribute 
in this direction, this work presents a careful redetermination of the Na abundances of 
the Hyades giants, motivated by the recent measurement of their angular diameters. 
An average of [Na/Fe] — +0.30, in NLTE, has been found. This overabundance can be 
explained by hydrodynamical models with high initial rotation velocities. This result, 
and a trend of increasing Na with increasing stellar mass found in a previous work, 
suggests that there is no strong evidence of Na overabundances in red giants beyond 
those values expected by evolutionary models of stars with more than ~ 2 Mq. 

Key vifords: open clusters and associations: individual: Hyades - stars: abundances 
- stars: evolution - stars: fundamental parameters. 



1 INTRODUCTION 

In many stages of their evolution, low- and intermediate- 
mass stars show signs of mixing between material of the sur- 
face with material of the interior that has been processed by 
nuclear reactions ('Pinsonncault"l997l; ICharbonnel fc TalonI 
L2OO8; Smil]anic_et al. 2009, and references therein). 

The standard model of stellar evolution, where con- 
vection is the only mixing mechanism, does not ac- 
count for all the observational details. The introduction 
of non-standard physical processes, such as atomic dif- 
fusion, rotation-induced mixing, internal gravity waves, 
magnetic buo yancy and thermohaline mixing, is unavoid- 
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As a star leaves the main sequence towards the 
red giant branch (RGB), its convective envelope deepens, 
causing first a dilu t ion of lithium, beryl li um, and boron 
jLebre et al.l 1 19991: IPasguini et al.l |2004| : [ Smilianic et al. 
I2OIOI: ICanto Martins et al.ll201ll ')~and then the first dredge- 
up ( Ibenlll967l ) , when material affected by hydrogen burning 
is mixed to the surface. The dredge-up causes an increase 
of the abundance of nitrogen and a decrease of carbon and 
of th e ^^C/^^C ratio (|Charbonnel et al.|[l99i : lGratton et all 
I2OOOI ). 

An extensive literature has shown a further mod- 
ification of the surface abundances after the bump in 
the luminosity function on the RGB. At this phase 
the abundances of Li, C and the ratio ^^C/^^C are 
further decreased while that of N increases. This ef- 
fect has been detected in stars of th e field and o f both 
open and glo b ular clusters (se e e.g. ISneden et al.l Il986l: 
Gilrovl Il989l: [GiTroy fc Brownl Il99ll: ICharbonnel et all 



1998 



ir iGratton et alj I2OOOI: iTautvaisiene et al.' '200GI. 
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1.1 The case of sodium 



2 THE HYADES 



During the first dredge-up, Na pr oduced by tiie NeNa-cycle 
ijPenisenltov fc Denisenl;ovalll990f ) in H-burning regions can 
potentially be mixed to the surface. The observational be- 
haviour of this element, however, is not clear. In open clus- 
ters, stars with Na overabundances as high as [Na/Fe|B 



-1-0.50 or more have been reported (iBraeaglia et al 



iFriel et al.ll2003l : IJacobson et al.|[2007i : ISchuler et al 



.20011 : 
2009^ ■ 



while others were found to have mild overabundances of 
[Na/Fel ^ -1-0.20 jHamdani et al.ll2000l: pTautvaisiene et al.i 



I2OO0I : iPasauini et al.ll2004l: iFriel et all 120101)" or abundances 
close to sola r ("Randich et al.' '2000; 'S estito et all l2007l : 
ISmihanic eTa l. 2009; Pancino ct al. 2010). 

One of the factors behind this discrepancy is the adop- 
tion of different log gfs for the same Na li nes by different 
autho rs. For example, the log gfs a dopted inlSmilianic et aP 
(|2009l '). from the NIST database (|Ralchenko et al.l I2OI0I) , 
are on average 0.22 dex higher than those adopted by 
Ijacobson et al.l (|2007l ). deri ved with respec t to A rcturus. 
Indeed, as discussed later bv lJacobson et all (|2008l ). a revi- 
sion of their gf-values resulted in an increase of about 0.20 
dex with corresponding decrease in the Na abundances. The 
use of this revised scale now results in modest overabun- 
danc es, [Na/Fe] = -h0.10...-h0.20 (see Table 14 of lFriel et af] 
l2010h in comparis o n wit h the strong overabundances found 
in lJacobson et all (|200it ). [Na/Fe] = -h0.40...-h0.60. 



Another issue are departures from the local thermody- 
namic equilibrium (LTE). Sever al authors provide non-LTE 
corrections for Na jBau muelle r et al. 199^ ; Gratton et al 



19991 ; iMashonkina et al.ll200CI;ITakeda et al.ll2003l ; IShi et al. 
2004l;fAndrievskv et al.ll2007l ; (Lind et al.ll201lh . As discussed 



bv lLind e t al.' bOlti. there is a scatter of 0.10 to 0.20 dex 
among the different corrections. Usually non-LTE abun- 
dances are smaller than the LTE ones. 

In the abundance analyses cited above, the different au- 
thors made different choices regarding these matters. It is 
thus difficult to make sense out of these results and under- 
stand whether there is a real problem with the sodium over- 
abundances in giants. Dedicated studies aiming to better 
understand this issue are still lacking in the literature. 

The re-analysis of Na in the Hyades giants presented 
here is a step on this direction. It is motivated by the recent 
determination of angular diameters (and thus fundamental 
effective temperatures - Toft) for the four Hyades giants 
from i nterferometric m easurements using the CHARA ar- 
ray bv lBovaiian et all (2009). Accurate Na abundances can 
help clarifying whether the strong overabundances are real 
and test whether there is an offset between observations and 
evolutionary models. 

This paper is divided as follows. Section [2] presents 
the observational data used in the analysis and the deter- 
mination of the atmospheric parameters of the Hyades gi- 
ants. Section [3] presents the determination of Na abundance 
from equivalent widths and spectrum synthesis while Sec. 
|4] presents a discussion of the results. Section [5] summarizes 
the conclusions. 



1 [A/B] = log [N(A)/N(B)]* - log [N(A)/N(B)]0 



The Hy ades are the closest open cluster to the Sun (~ 
46.5 pc, van Leeuwe n 2009). They have an age of 625 ± 
50 Myr JPerrvman et all 1 19981) and a metallicity of [Fe/H] 
= +0.13 ± 0.01 (|Paulsonet al.ll2003l 'l. There is no sign of 
Na overabundance i n the dwarf stars, i.e. [Na/Fe] ~ 0.00 
jPaulson et al.ll2003l ). For the g iants an average of [Na/Fe] 
= -1-0.48 (in LTE) was found bv lSchuler et all (|2009l '). which 
the authors regard as "too large to be explained by any 
known self-enrichment scenario" . 

The cluster has four evolved members, all He-burning 
clump giants; e Tan (HR 1409 or HD28305), 7 Tau (HR 1346 
or HD27371), 5^ Tau (HR 1373 or HD 27697), and 6^ Tau 
(HR 1411 or HD 28307). A fifth suggested gia nt member 
5 Ar i (HR 951) is likely a non-member ( do Bruiine et all 
l200j ). The stars e Tau and 7 Tau are s ingle stars while 5 
Tau and 9^ are s pectr oscopic binaries (|Griffin et al.l [l988l ; 
iMermilliod et all |2008| ). In addition, the star t Tau was 
found to have a massive planetary companion (|Sato et al.l 
[2007,) . 

An HR diagram of the Hyades on the region around 
the turn-o ff and th e clump is show in Fig. [T] An isochrone 
from Gira rdi et all (2002) with 625 Myr and [Fe/H] = 4-0.13 
indicates that a clump giant in the Hyades has ~ 2.48 M0. 
To illustrate, the corresponding isochrone is also shown in 
Fig. H] (no attempt to find the best fitting model was made). 



2.1 Observational data 

High-resolutio n spectra of th ree Hyades giants obtained with 
FEROS (Kaufer et al.lll999i ') at th e 2.2m MPG/ESO tele- 
scope at La Silla and with UVES jDekker et al.ll2000l ') fed 
by the UT2 of the VLT at Paranal are analyzed here. 

FEROS is a fiber-fed echelle spectrograph that provides 
a full wavelength coverage of AA350— 920 nm over 39 orders 
at R = 48 000. The spectra were reduced using the FEROS 
Data Reduction System (DRS) as implemented within ESO- 
MIDAS. UVES is a cross-dispersed echelle spectrograph able 
to obtain spectra from the atmospheric cut-off at 300 nm 
to ~ 1100 nm. Reducti on was done with th e ESO UVES 
pipeline within MIDAS jBallester et ahlbood ). 

The FEROS data of stars 6^ Tau and 7 Tau were 
made available to the author by Luca Pasquini (2010, pri- 
vate communication). The FEROS spectra of e Tau and the 
UVES spectra of e Tau and Tau were retrieved from the 
ESO/ST-ECF science archive facility. The log book of the 
observations is given in Table [1] 



2.2 Effective temperatures 

As the only giants of the nearest open cluster, these stars 
have been analyzed many times. To obtain an idea of 
the range of temperatures found in the literature, previ- 
ous determinations of the effective temperature (Tcff) of 
the four giants were queried at the PASTEL catalogue 
(|Soubiran et alj|2010t ). 

A few selected and recent results, together with a "lit- 
erature average" , are given in Table [2] along with the corre- 
sponding reference. N ote that the Teg calculated with the in- 
terferometric data bv lBovaiian et ah l (|20O9h is not included 
in the PASTEL catalogue. 
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Table 1. Observational data. 



Star 


Spectrograph 


V 


Date of 


Exp. Time 


S/N 








observation 


(s) 


@ 617 nm 


7 Tau 


FEROS 


3.654 


05 Mar. 2004 


120 


500 


e Tau 


FEROS 


3.540 


03 Oct. 2009 


180 


500 




UVES RED 580 




30 Nov. 2002 


2x1 


400 


5^ Tau 


FEROS 


3.764 


31 Oct. 2000 


120 


700 




UVES RED 580 




30 Nov. 2002 


2x1 


250 



Table 2. Effective temperatures of the Hyades giants taken from selected recent works from the literature. 



Star 


Tcfl 


Method 


Reference 




(K) 






7 Tau 


4844 ± 47 


Interferometry 


Bovaiian et al. (2009) 




4965 ± 75 


IRFMl 


Schuler et al. (2009) 




4960 ±8.1 


Line-depth ratio (average) 


Kovtvukh et al. (2006) 




4981 ± 80 


Average 


Selected results from the PASTEL catalogue 


51 Tau 


4826 ± 51 


Interferometry 


Bovaiian et al. f2009") 




4938 ± 75 


IRFM^ 


Schuler et al. (2009) 




4975 ± 7.6 


Line-depth ratio (average) 


Kovtvukh et al. (2006) 




5000 ± 80 


Fel excitation equilibrium 


Hekker & Melendez (2007) 




4968 ± 82 


Average 


Selected results from the PASTEL catalogue 


e Tau 


4827 ± 44 


Interferometry 


Bovaiian et al. f2009~) 




4911 ± 75 


IRFMl 


Schuler et al. (2009) 




4925 ± 8.7 


Line-depth ratio (average) 


Kovtvukh et al. (2006) 




4910 ± 80 


Fel excitation equilibrium 


Hekker & Melendez (2007) 




4925 ± 84 


Average 


Selected results from the PASTEL catalogue 



(1) The effective temperatures adopted bv lSchuler et al" I 1I2OO9I) were calculated bvlBlackwell fc Lvnas-Gravl lll994l) using the infrared 
flux method for 7 and 5^ Tau and bv lBlackwell fc Lvnas-Gra^Hl998t) for e Tau. 




Figure 1. HR diagram of the Hyades on the region around the 
turn-off and the clump. The four giants are shown with different 
symbols (e Tau as the blue square, 7 Tau as the green circle, <5^ 
Tau as the black triangle, and Tau as the red starred symbol). 
Luminosities and temperatures are from ide Bruiine et al.i (2001i) . 



There is a range of about 300-350 K o n the values of 
Teff d etermined for each star. As discussed in lBovaiian et all 
(|2009t ). the values determined from recent angular diameters 
are on the low side of this range. Similar ly, using the angular 
diameter of the Li-rich giant HD148 193. lBaines et al.l (|201ll ) 



derived an T^s that is on the cooler side of the range of 
temperatures determined in the literature. Together these 
results might be indicating that the temperature of giants 
is usually overestimated. 

In the following discussion, Na abundances are calcu- 
lated using two values of temp erature for eac h star: The 
interferometri c temperatures of Bovaiian et al.l (|2009l ) and 
those used bv lSchuler et al.l (|2009| ) in their abu ndance anal- 
ysis o f the same stars. The temperatures of iSchuler et al.l 
are close to the average literature values (see Table 
^ and differ from the interferometric ones by 90 to 150 K. 



2.3 Gravities 

Gravities for the Hyades giants were determined with the 
equation: 

log {9*/ge) = tog (M^/Mq) + 4log (Teff,, /Teff g) - log (L^/Lq) 



The luminosities bv lde Bruiine et al.l l|200ll ). a mass of 
2.48 M0, and the usual solar values (TbBq = 5777 K and 
loggQ = 4.44) were adopted. For consistency, gravities were 
calculated using each of the Teff- adopted for the analysis. 
In addition, to illustrate the effect of the gravity in the Na 
abundance, additional values were determined for 7 Tau as- 
suming masses of 2.0 and 3.0 M©. It can be seen that this 
change in the mass causes only a minor change on log g, ar- 
guing that this parameter is well constrained for these stars. 
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Table 3. The different sets of atmospheric parameters calculated 
for the Hyades giants. 



Star 


Tcff 


log g 


[Fe/H] 


Note 


7 Tau 


4844 


2.66 


+0.14 ± 0.05 


Interferometric T^g 




4844 


2.57 


+0.13 ± 0.05 


As above with 2.0 M 




4844 


2.74 


+0.15 ± 0.05 


As above with 3.0 M 




4965 


2.70 


+0.23 ± 0.05 


IRFM Tefi 


(51 Tau 


4826 


2.69 


+0.18 ± 0.06 


Interferometric T^ff 




4938 


2.73 


+0.25 ± 0.06 


IRFM Tcff 


e Tau 


4827 


2.60 


+0.26 ± 0.08 


Interferometric T^ff 




4911 


2.63 


+0.31 ± 0.08 


IRFM Tefi 



Table 4. Microturbulence velocities (in km s"^) derived with 
different calibrations for 7 Tau, using both the interferometric 
and the IRFM temperature and the corresponding log g. 



Calibration Interf. Teff IRFM Tpg 



© 


Gratton et al. (1996) 


1.37 


1.35 


© 


AUende Prieto et al. f2004) 


1.24 


1.35 




Carretta et al. (20041 


1.15 


1.15 




Alves-Brito et al. (2010) - T^tr 


1.26 


1.25 




Alves-Brito et al. ('20101 - loe e 


1.30 


1.29 



3 SODIUM ABUNDANCES 



2.4 Microturbulence 

Throughout the analysis, a fixed value of microturbulence, 
^ — 1.30 km s~^ was always ad opted. This is a t ypica l value 
found for the giants analyzed in lSmilianic et al.ll|2009l ). This 
value was checked against empirical relations from the lit- 
erature that calibrate ^ as a function of log g and/or Tcft. 
Five calibrations from four references were investigated. It is 
interesting to note that some of these calibrations are given 
only as a function of log g: 



(i) (, = 2.22 - 0.322 log g (|Gratton et al.lll996l . derived 
from giants with [Fe/H] betw een -1.00 to 0.00), 

(ii) C = 1-5 - 0.13 log g jCarretta et al.1 12004| . derived 
from open cluster giants with [Fe/H] between —0.37 to 
+0.24), 

(iii) C = 1-645 + (3.854 x 10^* (T^ft - 6387)) + ( 
-0.64 (logg- 4.373)) + (-3.427 x 1 0"^ {T^s - 6387) (logg - 
4.373)) (,Alle nde Prieto et al-lliooi . derived from solar neigh- 
borhood stars with [Fe/H] between —0.50 to +0.50), 

(iv) C = 3.40 - 4.41 X 10 ~^ Tetr and ^ = 1.84 - 0.202 
log g l|Alves-Brito et al.ll2010l . derived from bulge, thin and 
thick disc giants with [Fe/H] between —1.50 to +0.50). 

The values obtained from these calibrations for 7 Tau 
are given in Table |4] for the parameters using the interfero- 
metric and I RFM temperatures. A ll these calibrations, apart 
from that of ICarretta et all (|2004l ) give values that are very 
close to the one adopted here. What is more important, they 
show that the variation in Tcfi within the range considered 
here does not result in a large change in microturbulence. 
This argues that the choice of keeping ^ constant does not 
introduce systematic effects in the current analysis. 



2.5 Metallicity 

To estimate the metallicity of the stars ([Fe/H]), equiva- 
lent widths of a set of 15 selected Fe I lines were measured. 
The line list, atomic data, and equivalent widths are given 
in Table [H The Ce broadening constants were taken from 
ICoelho et al. I (|2005l). The ad opted solar iron abundance is 
A(Fe) = 7.50 (|Grevesse fc Sa uval 1998). Using the interfer- 
ometric temperatures, the mean metallicity of the Hyades 
giants is found to be [Fe/H] = +0.19 ± 0.06. Using the 
IRFM temperatures, this value increases to [Fe/H] = +0.26 
± 0.04. 



3.1 Line selection and atomic data 

All the 32 strong Nal and Nal l lines with wavelengths be - 
tween 4000 and 8200 A listed bv lSansonetti fc Martini (|2005| ) 
were considered as possible features to be used in this anal- 
ysis. The profil es of these lin es were checked both in the 
UVE^ and the lKurucd (|2005l ) solar spectra and in the spec- 
trum of 7 Tau. All lines that were heavily blended or too 
strong for an abundance analysis were discarded (e.g. lines 
at A 8183 and A 8194 A) .The atomic data of the remaining 
Nal lines are listed in Table [S] The excitation potential and 
the log fl.f of the lines w ere adopted from the NIST database 
(|Ralchenko et al" 201G|). The Ce broa dening constants wer e 
adopted from lCoelho et al] |2o2) and lBarbuv et ah l l|2006l ). 
An assessment of each line is given below: 

(1) 5 148.838 A: In the solar spectra the line has ~ 14 
mA. It is slightly blended on the blue wing with a line of 
similar strength. In the 7 Tau spectrum the two lines are 
completely blended but their bottoms can be distinguished. 
Only possible to analyze with spectrum synthesis. 

(2) 5 682.633 A: It is clear in the solar spectra that the 
blue wing is blended with a weaker line. In the spectrum of 
7 Tau the blend can not be recognized, but the line is asym- 
metric. Equivalent widths would be affected by the blend 
and thus spectrum synthesis should be preferred. In the Sun 
the line has ~ 100 mA. 

(3) 5 688.193 fc 5 688.205 A: These lines are the fine 
structure components of the same transition. The feature 
is strong in the Sun (> 120 mA) and stronger in 7 Tau (> 
170 mA). Seems to be clean enough, but it is too strong for 
an analysis using equivalent widths. 

(4) 6 154.226 A: On the Sun there is a weak line close to 
the red wing at ~ 6154.43 A. On 7 Tau the same blend is 
apparent but less distinguishable. There is another weak line 
at ~ 6154.6 A. The line seems to deviate from a Gaussian 
profile towards a Voigt one already in the Sun (where EW 
37 mA). Analysis using equivalent widths should be possible. 

(5) 6 160.747 A: Both in the Sun and in 7 Tau, the place- 
ment of the continuum is affected by the wing of the strong 
nearby Cal 6162.2 A line. It has closeby lines to the blue 
and red sides, but the profile seems clean. Analysis using 
equivalent widths is possible. It has ~ 60 mA in the Sun. 



^ The spectrum is available for download at the ESO website: 
www. eso.org/observing/dfo/quality/UVES/pipeline/solar .spec- 
trum. html 
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Table 5. Atomic data and equivalent widths of the Fe I lines used to derive the metallicity of the stars. 
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A X log gf Ce Sun 7 Tau Tau e Tau 

(A) (eV) (mA) (niA) (mA) (mA) 



5054.642 


3.64 


—2.032 


4.68E-32 


40.2 


80.6 


80.2 


84.5 


5127.679 


0.05 


-6.005 


7.38E-33 


22.3 


96.6 


99.5 


105.6 


5223.185 


3.64 


-2.285 


6.00E-32 


29.7 


65.8 


68.8 


76.0 


5320.035 


3.64 


-2.542 


8.91E-32 


19.8 


57.1 


57.9 


59.0 


5483.098 


4.15 


-1.481 


2.95E-31 


46.5 


78.8 


80.8 


83.1 


5522.446 


4.21 


-1.432 


3.02E-31 


43.6 


78.5 


80.6 


83.8 


5778.453 


2.59 


-3.524 


4.95E-32 


22.2 


73.6 


74.6 


81.1 


5784.658 


3.40 


-2.626 


3.57E-31 


27.2 


67.8 


70.7 


76.2 


5814.807 


4.28 


-1.861 


2.82E-31 


23.4 


56.5 


58.1 


59.7 


6012.209 


2.22 


-3.843 


3.35E-32 


24.2 


78.7 


78.4 


82.3 


6079.008 


4.65 


-1.055 


5.13E-31 


46.2 


77.8 


80.9 


81.2 


6187.989 


3.94 


-1.712 


4.90E-31 


47.9 


91.2 


93.1 


99.6 


6271.278 


3.33 


-2.797 


2.78E-31 


23.6 


68.3 


74.2 


79.5 


6739.522 


1.56 


-4.942 


2.10E-32 


11.8 


63.6 


64.8 


70.3 


6837.006 


4.59 


-1.756 


2.46E-32 


17.5 


47.1 


50.2 


52.1 



Table 6. Atomic data of the six selected Nal lines. 



A (A) 


log gf 


X (eV) 




5148.838 


-2.044 


2.102 


l.OlE-30 


5 682.633 


-0.706 


2.102 


3.38E-30 


5 688.193* 


-1.406 


2.104 


2.03E-30 


5 688.205* 


-0.452 


2.104 


3.38E-30 


6154.226 


-1.547 


2.102 


0.90E-31 


6160.747 


-1.246 


2.104 


0.30E-31 



(*) Fine structure components. 



3.2 Equivalent widths 

Given the above assessment, only the equivalent widths of 
the lines at 6154 and 6160 A were used to determine Na 
abundances. Howev er, other lines have b een used in the liter- 
ature. For example, ISchuler et all l|2009h also use the line at 
5682 A, which is clearly blended in the Sun and is asymmet- 
ric in 7 Tau. The equivalent width of such line should be re- 
garded as suspicious . Indeed, from the three lines adopted by 
ISchuler et alT(|2009l l, line 5682 A always results in an abun- 
dance that is higher by 0.10-0.14 dex than that obtained 
with the lines at 6154 and 6160 A. Without this line, the 
mean [Na/Fe] for the Hyades giants found by ISchuler et al.l 
l|2009l ) is reduced by 0.05 dex. 

Equivalent widths were determined by fitting Gaussian 
profiles to the observed lines using IRAE0. For the Sun, lines 
were measured both in the UVES and in the Kurucz spec- 
trum. The values obtained are listed in Table [T] 

Mo del atmospheres were computed usin g the Linux 
version l|Sbordone et all |2004| : ISbordond 120051 1 of the AT- 
LAS9 code originally developed by Kurucz (see e.g. iKurucj 
1 19931 ). For the calculations, the opacity distribution func- 
tions of lCastelh fc Kurucj (|2003l ') without overshooting were 
adopted. These models assume local thermodynamic equilib- 

IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



rium, plane-parallel geometry, and hydrostatic equilibrium. 
Abundances were derived using the WIDTH code, also in 
its Linux version. For the Hyades giants, abundances were 
calculated only with the equivalent widths measured in the 
FEROS spectra, as they have higher S/N. These Na abun- 
dances are listed in Table [S] 



3.3 Spectrum synthesis 

Abundances were also derived using spectrum synthesis and 
all the Na lines in Table |6l Synt hetic spectra were c alcu- 
lated with the codes described in ICoelho et al.l (|2005l ') and 
the model atmospheres described above. The line list is the 
one u sed to compute the spectrum library of ICoelho et al.l 
(120051 ). As mentioned before, the Na line at 5682 A is 
blended in its blue wing. In the line list used here, this blend 
is due to a Cr I line at 5682.495 A, with log gf = -0.609. 
It was modeled with the solar abundance recommended by 
iGrevesse fc Sauvall (|l998l ). A(Cr)© = 5.67. The Na abun- 
dances obtained with spectrum synthesis are listed in Table 

m 



3.4 Uncertainties of the Na abundances 

The main source of error of the abundances are the errors 
of the atm ospheric parameters. T he uncertainty of the Teff 
derived bv lBovaiian et ahl (|2009l ) is of the order of ± 50 K. 
This corresponds to an average uncertainty of ± 0.04 dex in 
the Na abundance - A(Na) - and ± 0.03 in [Fe/H]. 

As shown for 7 Tau (Table [21), gravities are well con- 
strained and no significant impact on the Na abundance is 
expected. A change of ± 0.05 dex in log g results in a change 
of ^ 0.015 dex in A(Na). The same uncertainty causes an 
effect of T 0.005 in [Fe/H]. 

The discussion in Section [2.41 shows that the value of ^ 
is well constrained. The values given by the different cali- 
brations have an rms of ± 0.08 km s~^. A total change of ± 
0.10 km s~^ in ^ results in a change of ^ 0.03 dex in A(Na). 
This change causes no significant effect in [Fe/H]. 

Considering that there is no significant uncertainty in 
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Table 7. Equivalent widths of three Nal lines measured from the data collected in this work. 



Line 


Sun 


Sun 


7 Tau 


51 Tau 


51 Tau 


e Tau 


e Tau 




Kurucz 


UVES 


FEROS 


FEROS 


UVES 


FEROS 


UVES 


(A) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


6154 


37.9 


37.7 


105.4 


104.1 


103.1 


108.0 


107.4 


6160 


58.1 


58.6 


119.6 


122.8 


119.8 


127.0 


124.8 



Table 8. Sodium abundances calculated with equivalent widths and the different atmospheric parameters of each star. 



Line 


Sun 


Sun 


7 Tau 


7 Tau 


7 Tau 


7 Tau 


<5i Tau 


51 Tau 


e Tau 


e Tau 


(A) 


Kurucz 


UVES 


Interf. 


IRFM 


2.0 M0 


3.0 M0 


Interf. 


IRFM 


Interf. 


IRFM 


6154 


6.27 


6.27 


6.86 


6.96 


6.86 


6.86 


6.83 


6.92 


6.91 


6.98 


6160 


6.32 


6.33 


6.81 


6.92 


6.81 


6.81 


6.86 


6.96 


6.95 


7.02 


Average 


6.29 


6.30 


6.84 


6.94 


6.84 


6.84 


6.85 


6.94 


6.93 


7.00 


[Na/Fe] 






+0.40 


+0.41 


+0.41 


+0.39 


+0.37 


+0.39 


+0.37 


+0.39 



the Na and Fe abundances of the Sun, the final uncertainty 
in [Na/Fe] is of ± 0.04 dex. 



4 DISCUSSION 

The mean solar abundance using equivalent widths is 
A(Na)0 — 6.30 while using spectrum synthesis it is A(Na)0 
— 6.33. The average [Na/Fe] values for each star are given 
in Tables [S] and [51 respectively for the cases using equivalent 
widths and spectrum synthesis. Although there is some dif- 
ference in A(Na) between the different temperature scales, 
there is basically no noticeable effect on [Na/Fe]. 

However, at first glance the average values determined 
with equivalent widths seem to be about ~0.12 dex higher 
than the values determined using spectrum synthesis. The 
difference is apparent in the giants but not in the Sun. 

As discussed below, this is not directly related to the 
use equivalent widths or spectrum synthesis to derive the 
abundances, but to the different choice of lines used in each 
case. There is however, a systematic difference of about ~ 
0.10 dex between the abundances derived with the A 6154 
line. This likely comes from an uncertainty in the equivalent 
width of this line. An error on the equivalent width of ± 
6mA can produce a change of ± 0.10 on the abundance. 

There are three effects contributing to the difference. 
One, is the Na abundance given by the A 5148 line. Exclud- 
ing it from the average, the [Na/Fe] values increase by ~ 
0.05 dex. It is not clear why this line results systematically 
in smaller values. Nevertheless, it was decided to consider 
its abundance as suspicious and to exclude it from further 
discussion. The final average LTE [Na/Fe] values, without 
line A 5148, are given in Table [TOl A second effect discussed 
in the next subsection are the NLTE corrections. 

The last effect seems related to differences between the 
WIDTH code, used to calculate abundances from equivalent 
widths, and the PFANT code, used to compute the synthetic 
spectrum. 

Following a suggestion by the referee, a series of syn- 
thetic spectra with the parameters of 7 Tau at the region 
around the line A 6160 were calculated with different Na 
abundances. The equivalent widths were measured and the 



values used to recompute the abundances. The resulting Na 
abundances given by WIDTH are 0.03 to 0.04 dex higher 
than the values used to compute the synthetic spectra. 

At this point, it is not possible to say whether this is 
caused by some numerical effect, by some difference in other 
input data (e.g. partition functions, opacities...), or to some 
difference in the physics, like the treatment of broadening, 
for example. This should, of course, be further investigated. 
Nevertheless, as it is based in the analysis of more features 
and likely based in a more reliable way to deal with the 
effects of broadening, the abundances derived with spectrum 
synthesis are preferred here. 

4.1 NLTE abundances 

As mentioned in the intr oduction, the N a lines are af- 
fected by NLTE effects (see lLind etallbOlll . and references 
therein). This has been suggested as a likely reason behind 
the larg e Na overabundances found in giants of op en clusters 
(see e.g. iRandich et al. I I2OO6I : ISestito et al]|2007l ). 

To correc t the abundances calculated here, the NLTE 
calculations of iLind et al.l (|201ll ) were adopted. Corrections 
were interpolated among the grid calculated in that work 
with an IDL routine kindly made available to the author 
by Karin Lind (2011, private communication). NLTE cor- 
rections were derived in a line-by-line basis, giving as input 
the atmospheric parameters and the LTE abundancefl 

Individual line corrections for all stars (giants and Sun) 
range from -0.03 (for line 5148 A) up to -0.15 (for line 6160 
A). For the Sun, the average corrections are of —0.08 and 
—0.09 dex, respectively with the abundances using equiva- 
lent widths and spectrum synthesis. Thus, in NLTE, the so- 
lar Na abundances derived in this work are A(Na)0 — 6.22 

* Alternatively, the interpolation code can accept as input the 
equivalent width of the line and return both the LTE and the 
NLTE abundances. The first appro ach was preferred here because 
the LTE abundances calculated in lLind et al.l 1I2OIII ) for a given 
set of atmospheric parameters are slightly different from the ones 
derived in this work. The difference is caused by different choices 
made in the treatment of line broadening. 
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Table 9. Sodium abundances with spectrum synthesis and the different atmospheric parameters of each star. 
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Line 


Sun 


Sun 


7 Tau 


7 Tau 


<5i Tau 


51 Tau 


e Tau 


e Tau 


(A) 


Kurucz 


UVES 


Interf. 


IRFM 


Interf. 


IRFM 


Interf. 


IRFM 


5148 


6.24 


6.24 


6.48 


6.56 


6.49 


6.58 


6.53 


6.62 


5682 


6.36 


6.36 


6.87 


6.97 


6.82 


6.92 


6.94 


7.00 


5688 


6.42 


6.42 


6.88 


6.98 


6.83 


6.94 


6.98 


7.05 


6154 


6.30 


6.30 


6.75 


6.86 


6.73 


6.83 


6.80 


6.88 


6160 


6.32 


6.32 


6.81 


6.91 


6.85 


6.96 


6.93 


7.00 


Average 


6.33 


6.33 


6.76 


6.86 


6.74 


6.85 


6.84 


6.91 


[Na/Fe] 






+0.29 


+0.30 


+0.23 


+0.27 


+0.25 


+0.27 



Table 10. Average [Na/Fe], in LTE, for each star and for each of 
the different analysis and parameters (excluding line A 5148 for 
the synthesis values). 



Table 11. Average [Na/Fe] in NLTE. 



Analysis 


7 Tau 


51 Tau 


e Tau 


EqW & Interf. 


+0.40 


+0.37 


+0.37 


EqW & IRFM 


+0.41 


+0.39 


+0.39 


Synthesis & Interf 


+0.34 


+0.28 


+0.30 


Synthesis & IRFM 


+0.35 


+0.31 


+0.32 



with equivalent widths and A(Na)Q =6.24 with spectrum 
synthesis. For the giants, the diflerent selection of lines re- 
sults on average corrections of —0.14 when using equivalent 
widths and around —0.11 when using spectrum synthesis. 

This difference in the average the NLTE corrections 
is another responsible for the difference among the LTE 
[Na/Fe] values given by equivalent widths and spectrum syn- 
thesis. With respect to the Sun, the correction when using 
equivalent widths (only lines 6154 and 6160 A) is of ~ —0.06 
dex. When using spectrum synthesis (4 lines), the correction 
is of —0.01 or —0.02 dex. This helps to explain why the dif- 
ferent NLTE abundances in Table[Tl]are in better agreement 
than the LTE values in Table [TOl 

The average absolute NLTE Na abundance of the 
Hyades giants, using the interferometric temperatures, is 
fo und to be [Na / Fe] = +0.3(0. With the IRFM adopted 
bv lSchuleret Zl (gooi) the value found here is [Na/Fe] — 

Using equivalent widths, ISchuler et al.l l|2009l) obtained 
[Na/Fe] = +0.48 in LTE. The question then is why the 
results are di f ferent . First, as pointed out in Section 13.21 
ISchuler et all l|2009l ) used the equivalent width of line A 
5682 to determine the Na abundance. This line, however, is 
blended and removing it from the average reduces [Na/Fe] 
by 0.05 dex. Second are the NLTE corrections for lines 
A 6154 an d 6160, causin g anot her reduction by 0.05/0.06 
dex. Last, ISchuler et all (|2009l ) adopted [Fe/H] = +0.13 
fo r the Hyades giants. The atmospheric parameters used 
m ISchuler et all |2009l) were determined m a previous pa- 



per. 



Schuler et al 



(I2OO6II. were Fel lines were also measured. 

Using these lines, ISchuler et al.1 l|2006l ) found an average of 
[Fe/H] = +0.16, although recalculating it with the Kurucz 
model atmospheres and codes used here a value of [Fe/H] 



Analysis 


7 Tau 


5^ Tau 


t Tau 


EqW & Interf. 


+0.34 


+0.33 


+0.32 


EqW & IRFM 


+0.35 


+0.33 


+0.34 


Synthesis & Interf 


+0.33 


+0.27 


+0.29 


Synthesis & IRFM 


+0.33 


+0.29 


+0.31 



= +0.19 is found. Taking into account this differe nce the 
final [Na/Fe] in NLTE found bvlSchuler et all ||2009|) should 
be [Na/Fe] = +0.34 (in the lSchuler et al.ll2006l ."metaUicity 
scale) or [Na/Fe] = +0.31 (in the metallicity scale recalcu- 
lated here). These values are in perfect agreement with the 
ones derived in this work. 



^ Using equivalent widths the value is [Na/Fc] 
^ Using equivalent widths the value is [Na/Fe] 



+0.33. 
+0.34. 



4.2 Comparison with evolutionary models 

An interesting question to look at now is whether the 
Na overabundances of the Hyades giants can be explained 
by evolutionary models. In standard models, no modifi- 
cation of the Na abundance is expected after the first 
dredge-up for stars below ~ 2.0 M© (|Mowlavil Il999l : 
ICharbonnel fc Lagardell2010l ). For stars of higher mass, an 
increase of up to +0.20 dex in the Na abundance is expected. 

When mixing induced by rotation is included in the 
models (transport of chemicals and angular momentum by 
shear turbulence and meridional circulation), larger over- 
abundances are produced. These effects can also create a dis- 
persion in the Na abundance among otherwise similar stars 
if they had diff erent initial rotation velocities in the zero age 
main sequence l|Charbonnel fc LagarddboiOl ) . This happens 
because rotation affects the internal abundance profile of the 
elements involved in H-burning. In this way, the Na-rich re- 
gion in rotating stars begins further out from the core, and 
more Na-rich material can be dredged-up to the surface. 

In Fig. [2] the Na abundance of the Hyades is 
shown in compar ison with the models calculated by 
ICharbonnel "fcLagardd (|2010i ). Also shown in the figure are 
the Na abundances of 31 giants of 10 open clusters derived in 
ISmilianic et~aLl |200i). These clusters have turn-off masses 
between 1.7 and 3.1 M0. The recommended Na abundance 
of the Hyades ([Na/Fe] — +0.30) is in the upper part of 
the range expected by the models. Therefore, and contrary 
to the conclusion of ISchuler et al.l (|2009l ), this comparison 
shows that the Na abundance in the Hyades can be explained 
by modern hydrodynamical models that include the effects 
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Figure 2. Sodium abundances, [Na/Fe], as function of tlie stellar 
mass at the clump. The circles in dicate the open cluster giants 
analyzed in lSmilianic et al] ||2009| ). the typical error bar of that 
work is shown in the lower right corner. The blue square corre- 
sponds to the recommended [Na/Fe] of the Hyades derived in this 
work. The lines represent the predicted [ Na/Fe] as a function of 
initial stellar mass given by the models of lCharbonnel &: Lagard3 
for the standard case (solid line), for a model including 
thermohaline mixing only (blue dotted line), and models with 
thermohaline mixing and rotation-induced mixing with initial ve- 
locities in the ZAMS of 250 and 300 km s~^ (lower red dashed 
line and upper green dashed line, respectively). 



of rotation. There is no need for an extra unknown mixing 
process. 

Although it is beheved here that the absolute Na abun- 
dance of the Hyades was derived, this claim can not be ex- 
tended to most of the results in the literature. One way to 
avoid such systematics and test whether the Na overabun- 
dances in giants conform with the prediction of the models, 
is to conduct a large homogeneous analysis of a sample in- 
cluding only giants but with different masses. In this way 
one can look for abundance trends with mass and test if 
they agree with the expectations of models. 

This is exempl ified by the giants analyzed in 
ISmilianic erahl ||2009|) As se en in Fig. [2] and discussed 
also i n Smilianic et al.l |200^ and ICharbonnel fc Lagardd 
there is an off-set of about 0.10 dex between obser- 
vations and models. However, there is an agreement in the 
increasing trend with mass. This suggests that we are indeed 
observing the effects of mixing in these stars, in spite of a 
possible systematic effect in the abundance scale. 

The Na abundance of the Hyades derived here and 
the trend of increasin g Na with increasing mass found in 
ISmili anic et al. I (|2009l ) argue that there is no strong evidence 
for overabundances above those expected by the models, for 
stars above 2 M©. 

At this point , it is interesting to mention the results 



of iPasauini et al.l ||2004| ) for the cluster IC 4651 and of 



iRandich et al. I l|2006l ) for M67. iPasauini et al] l|2004l ) found 



are not. Clump giants in M67 have ~ 1.3 Mq while clump 
giants in IC 4651 have ~ 1 .8 Mq. As can be seen in Fig. [2j 
according to the models of lCharbonnel fc Lagarda (120101 ) a 
star of 1.3 Mq is never expected to enrich itself in Na after 
the first dredge-up while stars of 1.8 Mq could be enriched 
by ~ -1-0.15 dex. Although caution is needed in the compar- 
ison among dwarfs and giants, these results seem to support 
the idea that models with rotation can properly explain the 
behaviour of Na also in giants below 2 Mq. 



4.3 Cluster versus field giants 

It is sometimes noticed that the overabundance of Na seen 
in gia nts of open clusters is not apparent in field stars l|Frie]| 
l2006l) . This comparison is, however, usuall y made between 
giant s in clusters and field dwarfs (see e.g. Ijacobson et al.l 
l201lh . As also discussed by these authors, such an offset 
might be real and caused by mixing. Based on the discus- 
sion of the previous S ection , this is a conclusion supported 
here. IJacobson et al.l (|201ll ) noticed indeed that their clus- 
ter abundances agree well with the Na abundances of field 
clump giants determined bv lMishenina et al.l l|2006l ). 

Comparisons among dwarfs and giants need to be care- 
ful not just because of mixing, but also because systematic 
effects can cause biases that might b e mistaken by real dif - 
fere nces. As an example, acco rding to lMelendez et al] (|2008l ) 
and lAlves-Brito et al] ||20ld) this seems to be the case be- 
hind previous claims of abundance differences among thick 
disc and bulge stars. 

A mismatch of Na abundan ces between st ar s in o pen 
clusters and in the field lead Ide Silva et al] (|2009l ) to 
suggest that the dissolution of open clusters might not 
be the main contributor of stars for the Galactic disc. 
These authors compiled Na abundances in cluster gi- 
ants from the literature, normalized them to a com- 
mon solar scale, and did a comparison wit h Na abun- 
dances in field dwa rfs from ISoubiran fc GirardI (|2005l ). field 
clump gi ants from Mishenina et al.l ||2006|). a nd bulge gi- 
ants from iFulbright. McWiUiam fc RichI l|2007l) . Agreement 
was found between the cluster giants and the field giants of 
J4ishcnina et al ( 20061). but an offset ex ists with respect to 
the dwarfs of ISoubiran fc GirardI (|2005l ) and the giants of 
IFulbright et all (|2007l ). 

Because of mixing, an offset between Na in giants and 
dwarfs might be expected, with the caveat that in samples of 
field stars one is adding together stars with different masses 
and metallicities, properties that affect the mixing of Na 
during the first dredge-up. In other words, depending on 
the mass range of the giants an offset between the Na abun- 
dances between giants and dwarfs and among giants them- 
selves might be expected or not. A r obust way to a t tempt 
a comparison such as the one done by Ide Silva et al] (|2009l ) 
would be using dwarfs in the field and dwarfs in clusters, 
where mixing is not able to affect the Na abundances. 



systematic difference at the level of 0.20 dex between the Na 
abundances of dwarfs and giants in IC 4651. A similar differ- 
ence between dwa rfs and giants was, ho wever, not detected 
in stars of M67 bv lRandich et al.l (|2006l ). Although these re- 
sults could be regarded as contradictory at first sight, they 



5 SUMMARY 

Sodium abundances of three Hyades giants have been rede- 
termined and an average value of [Na/Fe] — -1-0.30 in NLTE 
was found. This Na abundance was derived using the ab- 
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solute Tcfi of the stars determined with the interferomet ric 
angular diameter measurements of iBovaiian et all (|2009l ). 

This Na abundance agrees well with the ones predicted 
by th e hydrodynamical models of ICharbonnel fc Lagardd 
l|201Cl ) for a star of 2.48 M©, after the first dredge-up, and 
taking into account rot ation-induced mixing . This contra- 
dicts the conclusion of ISchuler et al.l l|2009l ) that the Na 
overabundances of the Hyades could not be explained by 
any known mixing mechanism. The Na abundance of the 
Hyades giants are on the upper-limit of the range predicted 
by the models, implying that the stars had a rather high 
initial rotation. 

Absolute abundance values have always to be consid- 
ered with care. Nevertheless, as a fundamental temperature 
was used, the Na abundance derived here should be quite 
accurate. In general, relative comparisons should be more 
robust. In this sense, the increasing trend of the Na abun- 
dance with incr e asing mass found in the giants analyzed by 
ISmilianic et~all (|2009i l is the same as the one expected by 
evolutionary mo dels. Agreement with models is also seen 
in the results of iRandich et al] (|2006f ). that did not find a 
difference in Na among dwa r fs an d giants of M67 and in 
the results of iPasguini et all (|2004l '). that did find a differ- 
ence in Na among dwarfs and giants of IC 4651. In addi- 
tion, it should be noticed that when comparing field giants 
with cluster g iants, similar Na overabundances are apparent 
l|Mishenina et a l. 2006.) . 

All these pieces of evidence seem to point to the conclu- 
sion that, so far, there seems to be no strong evidence for Na 
overabundances in giants of open clusters beyond those that 
can be well explained by the effects of evolutionary mixing, 
in stars more massive than ~ 2.0 Mq. A consistent and ho- 
mogeneous reanalysis of Na abundances in a large sample of 
giants is still necessary to confirm (or refute) this conclusion. 
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